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Abstract 

A A-invariant measure of a sub-Markov chain is a left eigenvector of its transition matrix 
of eigenvalue A. In this article, we give an explicit integral representation of the A-invariant 
measures of subcritical Bienayme-Galton-Watson processes killed upon extinction, i.e. upon 
hitting the origin. In particular, this characterizes all quasi-stationary distributions of these 
processes. Our formula extends the Kesten-Spitzer formula for the (l-)invariant measures of 
such a process and can be interpreted as the identification of its minimal A-Martin entrance 
boundary for all A. In the particular case of quasi-stationary distributions, we also present an 
equivalent characterization in terms of semi-stable subordinators. 

Unlike Kesten and Spitzer’s arguments, our proofs are elementary and do not rely on Martin 
boundary theory. 
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1 Results 


Let Z = (Z n ) n ^o be a subcritical Bienayme-Galton-Watson (BGW) process with offspring distri¬ 
bution of mean me (0,1). Denote by P the restriction of its transition matrix to N* = {1, 2,...}. 
Then P is a sub-stochastic matrix, the transition matrix of the sub-Markov process {Z killed upon 
hitting 0}. A measured u on N* is called a A-invariant measure for Z if it is a left eigenvector of 
P of eigenvalue A, i.e. if 

vP = Aza (1) 

In terms of generating functions, if F{z) denotes the generating function of the offspring distribution 
and G(z) = Yjk=\ u {k)z k the generating function of the measure u, then, supposing that G(z) is 
finite for all \z\ < 1 (a fact which follows from Lemma [TT] below), (jT]) is equivalent to 

G(F(z)) - G(F(0)) = AG(z), |z| < 1. (2) 


For xeN* denote by P x the law of the process Z starting from Zq = x and by E x expectation 
with respect to P x . Furthermore, for a measure v on N*, write |]zz|| = z/(N*). The following limit, 
called the Yaglom limit , is known to exist mm (see also El P- 16 ]): 


^min 


lim Pi(Z n e • 

n—► oo 


6 P n 

Zn > 0) = lim 1 

n—>oo OiP' 1 


( 3 ) 
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throughout the article, all measures are assumed to be locally finite unless explicitly stated. 

2 As usual, we consider measures as row vectors and functions as column vectors. 
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where the limit holds in the weak topology of measures on N*. Furthermore, the probability 
measure v m in satisfies © with A = m, i.e. it is an m-invariant probability measure of the process. 
In particular (see also [18, Proposition 5]), 


P!(Z n+1 >0) ||dlP ri+1 | 


m, asn-> oo. 


Pi (Z n > 0) Utfi P n \\ 

We denote the generating function of the probability measure by 


(4) 


H(z) = t'min (n)z n , \z\ ^ 1. 


n= 1 


Our main theorem is the following, which identifies all A-invariant measures of the BGW process 
(Z n )n^0- 

Theorem 1. 1. There exist no non-trivial (i.e. ^ 0) X-invariant measures for Z with A < m. 

2. The only m-invariant measures of Z are multiples of the Yaglom limit v m i n . 

3. Let a e (—oo,l). A measure v on N* is an m a -invariant measure for Z if and only if its 
generating function G(z) = ^ v(n)z n satisfies 

G(z)=\ (eW*)- 1 )* -e~ x ) — A(dx), \z\ < 1, (5) 

Jo xa 

where A is a locally finite measure on (0, oo) satisfying A (A) = A (mA) for every Borel set A a 
(0, go). The measure A is uniquely determined from v. Moreover, for every such measure A, (|5|) 
defines the generating function of an m a -invariant measure for Z with radius of convergence 
at least 1. 


Remark 2. In the proof of Theorem [lj the measure x~ a A(dx) will be constructed as the vague limit 
(n —* oo) of the measures p n on (0, oo) defined by fx n (A) = m~ an i , (p~ 1 A), where p n = Pi(Z n > 0) 
and A a (0, oo) Borel. 

Remark 3. We will give an overview over the existing literature in Section T3.31 but mention already 
here that Formula J5]) was obtained, in a slightly different form, for a = 0 and F(z) = 1 — m(l — z) 
(the “pure death case”) by Kesten and Spitzer [23j (giving credit to H. Dinges for deriving it 
independently). It was later shown by Hoppe [11] that the case of general and even multitype 
offspring distributions (but still a = 0) can be reduced to the pure death case. One could adapt 
Hoppe’s arguments for a ¥= 0, but we do not show this here. 


Quasi-stationary distributions. A A-invariant probability measure v (i.e., ||zz|| = 1) is also called 
a quasi-stationary distribution (QSD) of the process Z with eigenvalue A. The following result easily 
follows from Theorem [1] 

Theorem 4. A X-invariant measure v of the process Z is finite if and only if A < 1. In particular, 
v is a QSD with eigenvalue X of the process Z if and only if either 

1. X = m and v = v m i n , or 

2. X = m a for some a e (0,1) and the generating function G(z) = ^ v(n)z n satisfies (|5|) . where 
A is a locally finite measure on (0, oo) satisfying 

(a) A(A) = A(mA) for every Borel set A c (0, oo) and 

(b) $*(1 — e~ x )x~ a A (dx) = 1. 
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Furthermore, the measure A in (f5l) is uniquely determined from v. 

Remark 5. If A (dx) = \. dx in the above theorem, then G(z ) = 1 — (1 — H(z)) a , see Remark [TOl 
below. These QSD were found by Seneta and Vere-Jones in their seminal paper on QSD of Markov 
chains on countably infinite state spaces [22j . Rubin and Vere-Jones 1201 showed later that these 
QSD are the only ones with regularly varying tails and furthermore, every distribution v on N* 
with a tail of the form v([x, go)) = x~ a L(x) for a slowly varying function L(x) is in the domain of 
attraction of the above QSD (see also [21] for an analogous result for 1-invariant measures). 

Remark 6. The fact that the Yaglom limit i/ m i n as defined in ([3]) is the QSD of smallest eigenvalue 
is a general fact |6j, p515]. Furthermore, the fact that it is the unique QSD with eigenvalue m is 
classic in our case mm- 

Remark 7. QSD of BGW processes (and Formula Q) appear in a recent article by Henard and the 
author on random trees invariant under Bernoulli edge contraction [9]. 


Continuous-time BGW processes. A-invariant measures can be defined analogously for a sub- 
critical continuous-time BGW process (Z t )t^o- Let L and (Pt)t^o be its associated infinitesimal 
generator and semigroup, respectively, restricted to N*. We say that a measure v on N* is a A- 
invariant measure of the process (ZQt^o if z'L = — Xu, or equivalently, if vPt = e~ xt v for every 
t ^ 0. In this case, u is also a invariant measure of the embedded chain (Z rn ) n ^> o, for every 
r > 0. The measure A from Theorem [T] then satisfies A(A) = A(rA) for every Borel set A and every 
r > 0, hence A (dx) = - dx. We therefore have the following corollary to Theorem HI which also 
follows (in the pure death case) from results for general birth-and-death chains [4]. 

Corollary 8. Let (ZQt^o be a subcritical continuous-time BGW process and let m > 0 such that 
for all t ^ 0, P\[Zt] = e~ mt . 


1. 

2 . 

3. 


There exist no non-trivial (i.e. ^ 0) X-invariant measures for (Z t )t^o with A > m. 


The only m-invariant measures of {ZQt ^o are multiples of the Yaglom limit z/ m j„. 


For every A < m, the X-invariant measures of (Z t )t^o are 
whose generating function is given by Q with a = A/m 
explicitly, G is the generating function of a X-invariant 
c ^ 0, such that 


G(z) = c x 


r i-(i -H( Z )r, 

< -log(l - H(z)), 


exactly the multiples of the measure 
and A (dx) = \dx if X < m. More 
measure if and only if there exists 

a > 0 
a = 0 


(1 — H{z)) a — 1, a < 0. 


In particular, the only QSD of the process (Z t )t^o with eigenvalue am, a e (0,1] is the probability 
measure with generating function 1 — (1 — H(z)) a . 

Remark 9. A similar phenomenon happens for continuous-state branching processes, see m 

Remark 10. The explicit formulae in Corollary [8] are obtained from the following well-known equality 
which we recall for convenience: 


Va 6 (0,1) Va < 1 : 



dx 

x a+l 


Jr(-a)(a“ - 1), a ^0 

) — log a, a = 0. 


( 6 ) 


An easy proof goes by noting that for each a e (0,1), both sides of the equation define analytic 
functions on the half-plane {Rea < 1} and agree on {Rea < 0} as can easily be checked by 
calculating the two Euler integrals. The case a = 0 is also a special case of Frullani’s integral. 
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A-invariant measures of the process which is not killed at the origin. Say that a measure 
v on N = {0,1,...} is a true A-invariant measure for Z if it is a A-invariant measure for the non- 
killed process. In other words, if Pq denotes the transition matrix of the BGW process Z on N, 
a measure v on N is a true A-invariant measure for Z if and only if vPq = Xn, or equivalently, if 
its generating function G satisfies G(F(z)) = A G(z) for every \z\ < 1. Of course, since 0 is an 
absorbing state for the process, there are no true A-invariant measures for A < 1, and for A = 1 the 
only true (l-)invariant measures are the multiples of <5o (see e.g. [2j p. 67]). However, for A > 1 the 
A-invariant measures from Theorem [T] all extend to true A-invariant measures. In fact, we have the 
following analogue of Theorem [TJ 

Theorem 11. 1. There exist no non-trivial (i.e. ^ 0) true X-invariant measures for Z with 

A < 1. 

2. The only true (l-)invariant measures for Z are multiples of ho- 

3. Let a < 0. A measure u on N is a true m a -invariant measure for Z if and only if its generating 
function G{z) = ^n(n)z n satisfies 

f 00 1 

G(z)=\ eW*)- 1 )*—A(dx), |*| <1, (7) 

Jo x 

where A is a locally finite measure on (0, oo) satisfying A(A) = A (mA) for every Borel set 
A cz (0, oo). The measure A is uniquely determined from v. Moreover, for every such measure 
A, © defines the generating function of a true m a -invariant measure for Z with radius of 
convergence at least 1. 

For a subcritical continuous-time BGW process, we can define true A-invariant measures ana- 
loguously as above. The analogue of Corollary [8] is then the following: 

Corollary 12. Let (Z t )t^o be a subcritical continuous-time BGW process and let m > 0 such that 
for all t ^ 0, Ei [Z t \ = e~ mt . 

1. There exist no non-trivial (i.e. ^ 0/ X-invariant measures for (Z t )t^o with A > 0. 

2. The only true (O-)invariant measures for {Zt)tst o are the multiples of 5 o- 

3. For X < 0, the true X-invariant measures for {Z t )t ^o are exactly the multiples of the one 
given by (ED with a = A/m and A(dx) = \dx, i.e. the measures with generating functions 
G(z) = c(l - H(z)) a , c ^ 0. 

Overview of the article. The remainder of the article is organized as follows: Theorems [lj [4] and 
|TT]are proven in Section[2l Section[3]is an extended discussion consisting of the following three parts. 
In Section 13.11 we interpret Theorem [T] in the light of Martin boundary theory. Section 13.21 gives a 
probabilistic interpretation of the QSD from Theorem 2] in terms of semi-stable subordinators. In 
Section ECU we review the existing literature on A-invariant measures of BGW processes. 

Notation. Throughout the article, a statement involving an undefined variable z is meant to hold 
(at least) for every z 6 (0,1). 

Acknowledgments. I thank Olivier Henard for an extremely fruitful collaboration |9j from which 
this article arose. I also thank Alano Ancona and Vadim Kaimanovich for useful discussions about 
Martin boundary theory. An anonymous referee has made several valuable suggestions improving 
the presentation of the article. 
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2 Proofs 


We start with three simple lemmas. 
Lemma 13. Let z,w e (0,1), z A w. 

(l-p(l-z)) x b- 


Then for p > 0 small enough and for all x 5= 0, 

(l _ „)b/p ^ e (w-i)x _ -x if w < z 
1P> < e ifw>z 


Proof. Let z,w e (0,1), z A w. If w < z, then 1 — p( 1 — z) = 1 + p(z — 1) ^ e^ -1 ^ for all small 
enough p > 0. Furthermore, 1 — p ^ e~ p for all p > 0. This implies the first inequality. 

If w > z, then 1 — p ^ e f z ~ w ^ 1 )v f or a ll small enough p > 0. Furthermore, 1 — p( 1 — z) ^ 
for all p > 0. Hence, for p > 0 small enough and x ^ 0. 

(1 - p(l - z)) x/p - (1 - p) x/p ^ e (2 - 1)a: - e ( z ~ w ~ 1 ) x = e (z-w)x( e (w-l)x _ e -X} ^ e (w-l)x _ e -x 


This shows the second inequality and thus finishes the proof of the lemma. 


□ 


Lemma 14. Let v be an m a -invariant measure of Z, with aei. Set M = m 1 . Then for every 
(3 < a, there exists C < oo, such that 

v(\M n , M n+l )) < cW* n , Vn e N. 

As a consequence, we have for every (3 < a, for some C < oo, for every x ^ 1, 

jv([x, oo)) ^ Cx^ 13 , f3 > 0 
"|z/([l, x]) ^ Cx~P, /3 < 0 

In particular, XlnGN* o{n)\z\ n < oo for every \z\ < 1. Moreover, v is finite if a > 0. 

Proof. Fix a < M < A and e > 0. Let v be as in the statement and recall the definition of the 
transition matrix P. From the branching property and the law of large numbers we get for large n, 

5 n P([n/A,n/a)) = P n (Zi e [n/A,n/a)) > 1 - e. 

This implies that for every x large enough and y ^ x, 

vP(\x/A,y/a )) ^ (1 - e)v([x, y)). 

Now let £ —» 0. The previous inequality together with (jT]) (with A = m a ) gives for every (3 < a, for 
every x large enough and y ^ x, 


v([x, y)) ^ m fi v([x/A, y/a )). (8) 

Now set b n = v{[M n , M n+1 )). Iterating (j8j) and choosing A and a close to M one readily shows 
that for every /3 < a and 6 > 0, there exists I\ e N such that for n ^ K, 

bx+n m Pn (b K + • • • + |5n| +1)• 

Elementary arguments yield the first statement of the lemma. The remaining statements follow. □ 

Lemma 15. Let f : (0, oo) —> (0, oo) be measurable and satisfying for some constant C 1A 1: 

Vn e Z: either f = 0 on [m n , m 71 ^ 1 ) or Vx,y e [m n , m n_1 ) : / (x)/f(y) e [C _1 ,C]. 

Furthermore, let A be a measure on (0, oo) satisfying A(A) = A (mA) for all Borel A c (0, oo) and 
A([m, 1)) A 0. Then 


where we set a b 


1 


rcc 

f(x)A(dx) 
J o 


A (I>M)) Jo 

=> C~ x b < a < Cb. 


~c 
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Proof. First note that the restriction of the measure A([m, 1 )) -1 A to [m, 1) can be written as the 
image of the measure (logm^ 1 )” 1 ^ on [m, 1 ) under a suitable map (p: first map the latter via its 
distribution function to Lebesgue measure on [0,1], then map this back to [m, 1) via the inverse 
of the distribution function of the measure A([m, 1)) _ 1 A. Then extend the map (p to a map on 
( 0 , oo) by 

p{x) = m n lp(m~ n x), x e \m n+1 ,m n ), n e Z. 

By the self-similarity of the measures A and '-^f. the measure A([m, 1)) _1 A is indeed the image of 
the measure (log on (0, oo) by the map tp. Furthermore, by construction the map <p maps 

every interval [m”,m n_1 ), n e Z, to itself. In particular, for all x > 0, either f(tp(x)) = f(x) = 0 
or f (ip(x)) / f (x) e by assumption. The lemma easily follows by the change of variables 

formula. □ 


Corollary 16. Let p be a non-zero measure on (0, go) satisfying for some ael, p(A) = m. a p(mA) 
for all Borel A c; (0, oo). Then for every z e (0,1), 


rO 0 

(f 

Jo 


(H(z)-l)x _ - 


e x ) p(dx) < oo 


a < 1 . 


Proof. Note that /x([m, 1)) A 0, otherwise we would have p = 0 by self-similarity. Define A [dx) = 
x a p(dx). Then A satisfies the hypothesis of Lemma fT5l Now let / 3,7 e R u {+oo} and set 

fp,~i{ x ) = ^ 1 (^ 1 ) + x>0. 

By Lemma [T5l applied to the function x >—► x~ a fp^(x), we have for some C > 1 (depending on A, 
m, a, (3 and 7 ), 

rO 0 rO 0 rCO 

fp a (x)n(dx)=\ x~ a fg^(x) A(dx) x~ a ~ l fg n {x) dx. 

Jo Jo Jo 

In particular, this shows that 

1 

f/ 3 rf (x) p,{dx) < go (3 > a and 7 > a, 


J 

Jo 


(9) 


with the obvious meaning if/3 = +coor 7 = +ao. 

Let z e (0,1), so that H(z) e (0,1). We finally consider the integral 


rCC 

( e W*)-i)* - e -x)^ dx ). 
Jo 


( 10 ) 


For large x, the integrand is smaller than any fixed polynomial, so that the integral always converges 
at 00 by dHI) applied with /3 = +go and some 7 > a. On the other hand, as x —> 0, the integrand 
is asymptotically equivalent to H(z)x. Equation © applied with f3 = 1 and 7 = +00 then implies 
that the integral in (HOD converges at the origin if and only if a < 1. These two facts prove the 
corollary. □ 


Proof of Theorem{J\ Although we could restrict ourselves to the pure death case, i.e. F{z ) = 1 — 
m( 1 — z) (see Section Ill . we prove the theorem immediately in its generality. 

We first introduce some notation. Let Y n denote a random variable with the law of Z n under 
Pi(- | Z n > 0). By (©, Y n converges in law to the Yaglom distribution ^ m ; n , in particular, H n (z ) = 
—> H(z) as n —> 00 . Note that the inverse H~ l exists on [0,1] and is continuous. We further 
define p n = Pi (Z n > 0) for n e N and note that p n+ i/p n —* m as n —»■ 00 by (dJ). 
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Now let v be an m"-invariant measure, aeR. Denote by G(z) = X!n=i v ( n ) zTl its generating 
function, which is finite and well-defined for \z\ < 1 by Lemma 1141 We will extend the notation P x 
and Ej. to the (possibly infinite) measure u by P^(-) = n(n)P„(-) and Ej,[-] = ]F n u(n )E„[-]- 

Dehne the random variable N n to be the number of individuals at time 0 which have a descendant 
at time n. Then N n > 0 iff Z n > 0. Furthermore, by the branching property, Z n is equal in law 
(under P^ for every k e N*) to + • • • + Yn Nn \ where the variables are iid copies of Y n and 
independent of N n . Hence, as n —> co, by the m a -stationarity of n, 

m-^E^l^o] = m~ an E u [H~ 1 (z) Zn lz n >o\ = E u [H~\z) z °] 

= G{H- 1 {z)) G{H- 1 (z)), as n —> oo. (11) 

Now note that under P&, N n is binomially distributed with parameters k and p n for every 
k E N*. In particular, 

E 1/ [z JVn liv„ > o] = E u [z N "-0 N "] = E„[(l - Pn (l -z)) z ° - (1 -Pn) Zo ]. (12) 

Defining for every n e N the measure p n by p n (A) = m~ an v{p~ l A) for Borel A cz (0, oo), we thus 
get by (fill) and (fl2|) . 

r 00 

G(H~ 1 (z)) = hrn ((1 -Pn(l - z)) X/Pn - (1 -p n ) X/Pn ) Pn(dx). (13) 

n^co J Q 

With the first inequality in Lemma [13l and the fact that p n —> 0 as n —> go, this gives 

r 00 

Vw e (0,1) : sup ( e ( w ~ 1 ) x — e _a; ) p n (dx ) < oo. (14) 

n Jo 

Using (ED with w = 1/2, say, gives that the sequence of measures jl n (dx) = xe X p n (dx ) is tight 
and therefore, by Prokhorov’s theorem, precompact in the space of finite measures on [0, oo) en¬ 
dowed with weak convergence. Let p be a subsequential limit and dehne the measure p(dx) = 
x~ 1 e x Jl^ 0}O r ) ' ) (dx), where /!(o,oo) is th e restriction of p to (0,oo). We claim that 

r 00 

G(z) = (e^W- 1 )^ - e ~ x ) p(dx) + /2(0)i?(z). (15) 

Jo 

Indeed, fix z e (0,1) and denote by g z ,n(x ) the integrand on the right-hand side of (fT3l) . Then 
the function x >—> g z ^ n (x)/(xe~ x ), continuously extended to [0,oo), converges uniformly on every 
compact subset of [0, oo) to the function g z defined by g z (x) = ( e zx — l)/x for x > 0 and g z ( 0) = z. 
Using (|14p with some w e (z, 1) together with the second inequality in Lemma [13j a truncation 
argument then shows that we can pass to the (subsequential) limit inside the integral in (1131) , which 
yields 

CO 

g z (x)p(dx). 

) 

This yields (1151) . Furthermore, the theory of Laplace transforms gives that p and g(0), hence p, 
are uniquely determined by (fT5j) . so that p n converges in fact weakly to p. As a consequence, p n 
converges vaguely on (0, oo) to p. 

The scaling properties of the measure p follow from this convergence: we have for every compact 
interval A cz (0, oo) whose endpoints are not atoms of p , 

p(A) = lim ?n _a ( n+1 (u(p“! 1 A) = m~ a lim m~ an v{jp~ 1 (p n /p n+ i)A) = m~ a p(m~ 1 A), 
n —>oo " r n—► oo 

since p n +i/p n rn by ©• This implies that the measure A defined by A (dx) = x a p(dx) satisfies 
A(A) = A (mA) for every Borel set A. 


G(H-\z)) = f 
Jr 


7 


It remains to investigate which terms in (11511 vanish for particular values of a. A first constraint 
comes from the fact that G(z ) is finite for every ze (0,1) by Lemma ITU and so the integral in (1151) 


needs to be finite as well. By Corollary 1161 this is true if and only if a < 1 or /u = 0. 

A second constraint comes from the fact that G satisfies ([ 2 ]) with A = m a . To verify this, we 
first recall the following equations for the function H: 

H(F(z)) - H(F(0)) = mH(z), \z\ < 1. (16) 

H(F(0)) = 1 — m (17) 

H(F(z)) - 1 = m(H{z ) - 1), \z\ ^ 1. (18) 

Indeed, (11611 is an immediate consequence of © (with A = m) and the finiteness of H for |z| = 1, 


(El) follows from (1161) by setting z = 1, and (1181) follows from © and ED by reordering terms. 
We now have by (fl5l) . for every 2 £ (0,1), 


G(F(z)) - G(F( 0)) 

rCO 

[ e (H(F(z))~l)x _ e -x _ e (H(F( 0))-l)x + e -xj + ^(Q)( F (~)) _ H (F( 0))) 

Jo 

rCO 

[ e m(H{z)-I)x _ e mx j + ( z ) ( by EJ , ED , ED (in this order)) 

Jo 

r oo 

1 [e(' ff ( 2! ) -1 ) x — e x ] fj,(dx ) + mfl(0)H(z) (by self-similarity of fi). 

Jo 


= m 


Comparing with (J2J) , this implies that p,( 0) = 0 unless a = 1. Summing up, we have the following 
constraints for the quantities in () 15 1) : 

• a > 1 : fi = 0 and p,(0) = 0 


• a = 1 : /i = 0 

• a < 1: /x(0) = 0. 

This proves the necessity part of the theorem. 

For the sufficiency, we only need to consider the case a < 1. Let G be a function given by 
© with A a measure on (0, oo) satisfying A(A) = A (rnA) for every Borel set A. By the above 
calculations, one readily shows that G satisfies © with A = m a . It remains to show that G is the 
generating function of a (locally finite) measure on N*. Now, for every x e (0, oo), the function 
z >—> e^G)- 1 ) x _ e ~ x i s generating function of the sum of M iid random variables distributed 
according to v m \ n , where J\f ~ Poi(x), restricted on the event that this sum is positive (see also 
Section T3. 2D . Hence, G is an integral over a family of generating functions and thus the generating 
function of a (not necessarily locally finite) measure. But by Corollary El G(z) is finite for z £ ( 0 , 1 ), 
so that this measure is indeed locally finite. This finishes the proof of the sufficiency part of the 
theorem. □ 


Proof of Theorem [7} Let v be a non-trivial m"-invariant measure of the BGW process Z, ct < 1. 
By Theorem [H it remains to show that v is finite if and only if a £ (0,1]. For a = 1 this is 
immediate, suppose therefore that a < 1. Denote by G the generating function of the measure 
v and let A be the measure from Theorem |T] Then Lemma El easily implies that the integral 
^°(1 — e~ x )x~ a A(dx) converges at the origin for all a < 1 but converges at co if and only if 
a £ (0,1]. Hence, ||^|| = G(l) < oo if and only if a £ (0,1]. This proves the theorem. □ 

Proof of Theorem lill The first two parts are known, see the discussion before the statement of the 
theorem. The third part can be proven by adapting the proof of Theorem [33 Alternatively, it can 


be derived from Theorem [l] as follows: Let A = m a > 1 (hence, a < 0). Let v be a measure on N 
and denote by v* its restriction to N*. Denote by G and G * the generating functions of v and v*, 
respectively, note that G* = G — G( 0). Since 0 is an absorbing state for the process Z, the measure 
v is a true A-invariant measure for Z if and only if the following two statements hold: 

1. v* is a A-invariant measure for Z. 

2. zaP(O) = Az'(O), equivalently, G(F( 0)) = AG(0). 

By Theorem [U the first statement is equivalent to 

f 00 1 

G(z)=\ (dx)+C, 

Jo x ° 

for some constant C and A as in the statement of Theorem |T] (it can easily be seen that the 
integral converges using Lemma fl5l as in the proof of Corollary fT 6 |) . Together with (|17j) and the 
self-similarity of A, this gives 

G(F( 0)) = f ” e~ mx — A (dx) + C = A f°° e~ x — A (dx) + C = AG(0) + C. 

Jo xDl Jo 

Hence, given the first statement, the second statement is equivalent to C = 0, which proves the 
theorem. □ 


3 Discussion 

3.1 The Kesten—Spitzer formula for invariant measures and the minimal Martin 
entrance boundary 

To our knowledge, Theorem |TJ and more specifically Formula (]5|) , was previously known only for 
A = 1 (i.e., a = 0). In this case, one simply says invariant instead of 1-invariant. In the literature 
(Kesten-Spitzer [23], Athreya-Ney [2] p. 69], Hoppe [11] ; see Section (T3] below for the history of the 
result), one generally finds this result under the following form: A function Q(z) is the generating 
function of an invariant measure for the BGW process Z if and only if there exists a constant c ^ 0 
and a probability measure /r on [ 0 , 1 ), such that 

rl 00 

Q(z) = c ^ [exp((Lf(;z) — l)m n_t ) - exp(-m n ^*)] n(dt). (19) 

JO n —_QQ 

This is the Choquet decomposition of Q(z ) as a convex combination of generating functions of 
extremal invariant measures. One easily sees that ([5]) (with a = 0) and (|19D are equivalent: Given 
c ^ 0 and a measure fj, such that (fT9|) holds, we can define a measure A on [l,m -1 ) as the push- 
forward of the measure c/i by the map t >—> m - L The measure A can then be uniquely extended to 
(0, oo) in such a way that A(A) = A (mA) for every Borel set A. One easily checks that ([5]) holds 
with this A and a = 0. Conversely, given such a measure A, one can define a finite measure p, on 
[0, 1) as the push-forward of the measure A(- n [1, m _1 )) by the inverse map x <—> log m -i x. Setting 
c = /i([0, 1)) and /r = p/c gives (fl9l) . 

We now relate formula (1191) to Martin boundary theory, see m Chapter 10] for an introduction 
to this theorjU. We briefly recall the basic constructions of interest to us. Let P be the transition 
matrix of a transient sub-Markov chain on N*. Define the Green kernel G(x,y) = Xi^=o (P k )xy and 

3 Another very good and more modern introduction is [24l Chapter IV]. He only considers Martin exit boundaries 
but one can reduce to this case in our setting by considering the transition matrix Pij = v / ln in(j)Pji/v I nin{i) instead of 
P. 
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assume there is a state o e N* such that G(x, o) > 0 for all x e N* (in the case of subcritical BGW 
processes killed at 0, we choose o to be the span of the reproduction law). This allows to define the 
Martin kernel by K(x,y) = G(x,y)/G(x,o). The Martin entrance compactification of N* is then 
defined as the smallest compactification A4 of the discrete set N* such that all measures K(x,-) 
extend continuously (w.r.t. pointwise convergence of measures seen as functions on N*). Every 
point £ on the Martin entrance boundary B = ,A/f\N* thus defines an invariant measure A"(£, •) with 
mass 1 at o. Moreover, every extremal invariant measure, meaning that it can not be written as 
a non-trivial convex combination of invariant measures, arises this way. The set of those points 
£ e B for which A(£, •) is extremal is called the minimal Martin entrance boundary, denoted by B*. 
The Poisson-Martin integral formula now assigns to every invariant measure v a unique integral 
representation in terms of extremal invariant measures, namely, 

V=\ A'(£,')/m), 

JB* 


for a finite measure pP on B*. 

The construction outlined in the previous paragraph is the approach used by Kesten and Spitzer 
[23] to derive formula (fl9l) (for pure death processes). In particular, their proof implies that the 
extremal invariant measures of a subcritical BGW process are (up to multiplicative constants) the 
measures u t , t e S 1 = [0, l]o~i, with generating functions 

00 00 

^ vt(k)z k = ^ [exp {{H(z) - 1 )m n ~ t ) - exp(-m n ~*)]. 

k =1 n =—oo 

Defining £ f e B * by v t = A"(£ f ,-), the map t >—> £t is thus (by extremality) a bijection between 
the compact space S 1 and B*, moreover, one easily sees that it is continuous. It follows that the 
minimal Martin entrance boundary B* is homeomorphic to the circle S 1 . 

Now let A > 0. The above construction can be performed with the operator P/A instead of P, 
giving rise to a A-boundary theory for all A such that the A-Green function G\ = X!*=o A ~ k {P k ) X y is 
finite. The infimum of these values of A is the spectral radius p = limfc- >00 ((P fc ) 00 ) 1//fc [Ml Chapter II], 
which equals m for subcritical BGW processes by Q and (JH). For A > p = m, Theorem [T] then 
implies a formula similar to f)19[i . A reasoning as in the last paragraph yields the following: 

Corollary 17. For every A > m, the minimal X-Martin entrance boundary of the BGW process 
(Z n ) n ^Q is homeomorphic to the circle S 1 . 

In particular, Corollary [T7] shows that all minimal A-Martin entrance boundaries, A > m, are 
homeomorphic. This remarkable fact is part of a property called stability by some authors [Ml p301] 
and holds true for example for the (exit) boundary of random walks on trees and hyperbolic graphs. 
We know of no general theory that yields this result without explicitly calculating the A-Martin 
entrance boundaries for every A. 

We finish this section with a discussion of the case A = m, for which Theorem [T] gives that 
the minimal m-Martin entrance boundary is trivial, i.e. there exists up to multiplicative constants 
only one m- invariant measure. This fact is quite common and holds in general for example if the 
process is m-recurrent, i.e. if G m (x,y) = oo for all (some) x,y [M, Chapter IV], Note that m- 
recurrence is equivalent to recurrence of the so-called Q-process, which is in our case the Markov 
process with transition matrix Q given by Qij = m~ 1 P{jj/i (recall that the function h(i) = i is 
m-harmonic for our process, i.e. Ph = mh). It is remarkable that in our setting the Q-process 
may be positive recurrent, null recurrent or transient. This fact does not seem to appear in the 
usually cited monographs on branching processed, only a criterion for positive recurrence is easy 

4 It was even claimed in the literature that recurrence always holds 1191 p972]. 
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to find (see e.g. |2j p59]): the Q-process is positive recurrent if and only if Ei[ZilogZi] < oo. 
However, Joffe proved in 1967 already the following recurrence criterion m- Let F(z) denote the 
generating function of the offspring distribution and define 77 by 1 — F(z) = m( 1 — z)( 1 — 17 (z)). Set 
q n = P 1 (Z n = 0). Then the Q-process is recurrent if and only if the following sum diverges: 

00 n 

2 na-*«))■ 

n= 1 k =1 

Since 1 — q n = m n+ °^ by (jU), one can easily construct examples where the above sum converges 
(so that the Q-process is transient), for example when rj(z) > 1/| log(l — z)\@ for some (3 < 1 and z 
close to 1. 

3.2 Probabilistic interpretation of (J5]) and relation with semi-stable subordina- 
tors 

Let v be a QSD of eigenvalue m a of the BGW process, a e (0,1). By Theorem 0] it admits the 
representation d2|) with a measure A as in the statement of the theorem. Let Af be a random variable 
whose generating function is equal to the right-hand side of ([5]), but with H{z) = z. Then v is the 
law of the sum of Af iid random variables distributed according to the Yaglom distribution z/ min . 
As for the law of Af, expanding the exponential in © gives 

\/k ^ 1 : P(Af = k) = [ e~ x ^-^A(dx), P(AA = 0) = 0. (20) 

Jo x a 

Heuristically, Af is therefore a Poisson-distributed random variable with a random parameter drawn 
according to the measure x~ a A (dx) and conditioned to be non-zero. A way to make this rigorous 
(note that the measure x~ a A(dx) has infinite mass!) is using subordinators, of which we first recall 
the basic facts. 

A subordinator S = (St)t^o is a real-valued, non-decreasing process with stationary and inde¬ 
pendent increments. We always assume So = 0. Then the law of S is determined by its cumulant 
ks{0) = — log Efe -05,1 ], which satisfies the Levy-Khintchine formula (see e.g. [HI Ch. 13] or [3]), 

POO 

k s ( 0) = a 6 + (1 - e~ ex )M{dx), (21) 

Jo 

where a ^ 0 is called the drift and M is a measure on (0, 00 ) called the Levy measure of the 
subordinator and satisfying $®(1 — e~ x )M(dx ) < 00 . 

If T = (Tt)t ^0 is another subordinator (or, in general, a Levy process) independent of {St)t^ o> 
then the subordinated process ToS := (Ts t )t ^0 is again a subordinator (Levy process) with cumulant 

KToS = K S 0 KT- ( 22 ) 

If N = (N t )t^o is a driftless subordinator whose Levy measure is a probability measure on N*, 
then the subordinators N and N o S both take values in N* and N o S is therefore again a driftless 
subordinator with Levy measure concentrated on N*. Let H and G denote the generating functions 
of the Levy measures of N and N o S, respectively. Note that 1) = 1 and G(l) < 00 , because a 
Levy measure on N* is necessarily finite. It follows from (12111 (applied first to N o S and then to 
N) and (|22l) that 

G(l) - G(z) = k N oS (- log z) = k s (1 - H{z)). (23) 

Setting z = 0 yields G(l) = /ts(l). Rearranging ([23]) . we get with (|2Tj) . 

POO 

G(z) = k s (1) ~ « s (l - H(z)) = aH(z) + ( e ^~ 1)x - e~ x ) M{dx). (24) 

Jo 
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We apply the previous equations to the QSD u, by setting a = 0 and M(dx ) = -^A(dx), 
note that ks( 1) = $*(1 — e~ x )M{dx ) = 1. In particular, M is a Levy measure, so that the 
subordinator S is well defined. We also let the Levy measure of the subordinator N be z/ m i n ; 
we recall that its generating function is indeed denoted by H(z). Equation (1231) then gives a 
probabilistic interpretation to the QSD v. it says that v is the Levy measure of the subordinator 
N o S (or, equivalently, the law of its first jump). 

Note that the case a = 1 may also be covered by setting a = 1 and M = 0 in (1211) . i.e. taking 
the subordinator S = Id. 

This fact allows for an alternative statement of Theorem [3J For this, we introduce the notion of 
a semi-stable subordinator: we say that the subordinator S = (St)t^o is ( a,m)-semi- stabl^, if 

o L = (mSt) t> o, (25) 

or, in terms of the cumulant, 

Ks(m9) = m a Ks{0 ), 9^0. (26) 

One easily obtains from ()21l) and (|26l) the following characterization of semi-stable subordinators: 
A subordinator S = (St)t^o with drift a, Levy measure M, satisfying Sq = 0 and S ^ 0, is 
(ct, m)-semi-stable, a e K, m e (0,1), if and only if 

• a e (0,1), a = 0 and M(dx) = A (dx) for a measure A on (0, oo) satisfying A(A) = A (mi) 

for all Borel A cz (0, oo), or 

• a = 1, a > 0 and M = 0. 

The previous arguments then give the following equivalent statement of Theorem [4j 

Theorem 18. The quasi-stationary distributions of eigenvalue m a of the BGW process, ael, are 
exactly the Levy measures of the subordinators N o S, where N is the driftless subordinator with 
Levy measure v m i n and S is an (a, m)-semi-stable subordinator with ks( 1) = 1. 

Remark 19. One can drop the requirement ks( 1) = 1 in the above theorem if one replaces “are 
exactly the Levy measures” by “are exactly the laws of the first jumps”. 

Composition of generating functions Let G a be the generating function of a QSD of Z with 
eigenvalue m a , a 6 (0,1]. Furthermore, let Gg be the generation function of an m“^-invariant 
measure, j3 ^ 1, of the pure death process with mean offspring m a , i.e. with F(z) = 1 — m“(l — z). 
It is easy to see from ([2]) that the composition GpoG a is the generating function of an m a P- invariant 
measure of Z (note that the Yaglorn distribution of a pure death process is always <5i, hence its 
generating function is the identity z >—► z). If A a , A^ and A a p are the measures from Theorem [1] 
corresponding to G a , Gp and Gp o G a , respectively, then one may ask the following question: 

Question 20. Is there a simple formula expressing A a p in terms of A a and Ap ? 

We were not able to answer this question and are in fact doubtful that the answer is positive in 
general. In order to rephrase this problem into a more familiar setting, consider the case where Gp 
is the generating function of a probability measure, so that in particular /3 6 (0,1]. Let S a and SP 
be the (a,m)- and (/3, m a )-semi-stable subordinators associated to G a and Gp by Theorem [T8l In 
particular, K 5 «(l) = /c^Q) = 1. By fl23l) and (l22l) . we then have 

1 -G^o G a (z) = K S p( 1 - G a (z )) = k s p(k s <x( 1 - H(z))) = n SaoSI s(l - H{z)). 

Hence, Question [20] is equivalent to the question of whether there is a simple formula expressing 
the Levy measure of S a o S 13 in terms of the Levy measures of S a and . To the best of our 
knowledge, no such formula is known, and, given the fact that the Laplace transform of a measure 
has no simple inversion formula, there does not seem to be much hope. 

5 This terminology is taken from [5} Section 9.2]. 
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3.3 History of the problem 

The study of A-invariant measures of subcritical BGW processes has a rich history which we aim 
to elucidate here. The starting point seems to be Yaglom’s 1947 article [25], who showed the 
existence of the now-called Yaglorn limit of a subcritical BGW process under the assumption of 
finite varianc^. The BGW process appeared again as an important example in the seminal paper 
by Seneta and Vere-Jones [22] on QSD of Markov processes on (countably) infinite state spaces. 
In this work, the authors show that subcritical BGW processes admit a one-paranreter family of 
QSD whose generating functions are 1 — (1 — H(z)) a , a e (0,1], with H(z ) denoting, as above, the 
generating function of the Yaglorn limit. Rubin and Vere-Jones m raised the question whether 
there existed other QSD. They failed to answer the question in general but showed that these QSD 
where the only ones with regularly varying tails. 

These works on QSD of subcritical BGW process seem to have been independent of other works 
on (l-)invariant measures: In 1965, Kingman [T6] showed that invariant measures for a subcritical 
BGW process are not unique, which, as claimed by Kingman, disproved a conjecture by Harris. A 
full characterization of invariant measures, Formula (1191) . was then given by Kesten and Spitzer in 
1967 [23] (they also gave credit to H. Dinges for deriving the formula independently), motivated 
by the need of finding examples of explicitly calculable Martin boundaries for Markov processes. 
Spitzer’s note only contained a brief sketch of a proof and covered only the pure death case, but he 
claimed that the method would work as well for arbitrary offspring distributions if E[Zi log Z i] < go. 
A full proof of this fact appeared in Athreya and Ney’s well-known monograph [2, p. 69], which also 
covers the Yaglorn limit but does not treat QSD in general. 

In the 1970’s, Hoppe considered again the question of the uniqueness of the QSD with generating 
functions 1 — (1 — H(z)) a , a e (0,1]. As many of the previous works on branching processes, he 
extensively used generating functions. Starting point was the following equation, which, for the 
generating function G of a probability measure v, is easily seen to be equivalent to (|2|) : 

l-G(F{z)) = X{l-G(z)). (27) 

Hence, finding all QSD of eigenvalue A amounts to finding all probability generating functions G 
solving ([271) . Hoppe [10] showed in 1976 that one can reduce the problerrQ to the pure death case 
F(z) = 1 — m( 1 — z): He proves that a generating function G satisfies (1271) with A = m a if and only 
if there exists a generating function A(z), such that G(z) = A(H(z)) and 

1 — A{ 1 — m( 1 — z)) = m a ( 1 — A(z)). (28) 

He also remarks that the general solution A(z) to this equation is of the form 

A(z) = 1 - (1 - z) a expty(- log(l - z))), (29) 

for a [ log m |-periodic function with tjj(0) = 0. The drawback of this representation, apart from its 
uncertain probabilistic meaning, is that it is not immediate from (1291) whether the Taylor series of 
A(z) only has non-negative coefficients, i.e. whether A(z) is the generating function of a probability 
distribution. Hoppe [IT5] was not even sure whether such a function exists for a non-constant i/j. 
However, one can show (using for example theorems by Flajolet and Odlyzko [7, Proposition 1[) 
that for every c ±,..., c n there exists Co > 0, such that for |c| < Co, the Taylor expansion at 0 of the 
function 

A(z) = 1 - (1 - z) a exp c k sin log(l - z)^j 

®The assumption of finite variance was later removed in M- 
7 He also showed in another article m that this is true for invariant measures as well, which allowed him to prove 
Formula m without additional conditions on the offspring distribution. Note that Formula (HU) was again reproven 
in the general case in |T], the authors of which were apparently unaware of Hoppe’s work. 
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only has non-negative coefficients (a similar reasoning has been used by Kingman in his article cited 
above [16)). The function is therefore a generating function of a probability distribution which is a 
QSD of the pure death process. This gives an alternative proof of non-uniqueness of the QSD but 
no satisfying characterization. 

In 1980, Hoppe [12] therefore published another representation of solutions of ([28]) : He showed 
that there exists a one-to-one correspondence between QSD and invariant measures of the BGW 
process. Again, he used functional equations: by ([2]), a (non-trivial) measure v on N is an invariant 
measure of the BGW process if and only if there exists a normalizing constant c > 0, such that the 
generating function Q(z) = Yjh=\ cn^i)z n satisfies the functional equation 


Q(F(z)) = 1 + Q(z), Q( 0) = 0. (30) 

Hoppe [12] then showed that for every a e (0,1], the functiord 

G a (z) = ¥ - — -——- (31) 

H' (w)m^ a ~ 1 lDM dw 


is the generating function of a QSD of eigenvalue m a of the BGW process and conversely, for every 
such function, setting 


Q(z) 


!og(l - Gg(z)) 
log m a 


(32) 


defines a generating function which solves (1301) (note that this is a special case of the compositions 
of generating functions studied at the end of Section 13.21) . This yields for every a e (0,1) a 
bijection between all QSD of eigenvalue m a and all invariant measures and thus apparently solves 
the problem of characterizing all QSD. However, the non-linear transformations from Equations 
(13TT) and (1321) do not seem to be easy to tame, for example, we are not aware of any direct way of 
obtaining a formula like ([5]) from (1191) using the above formulae. More specifically, we are unable 
to relate the measures A in the respective representations of G a and Q in ([5]), when G a and Q are 
related through (1311) or ([32]) . We do not believe that there exists a simple relation between them, 
similarly to our reservations concerning Question [20] Therefore, to the best of our knowledge, the 
current article provides a new approach to A-invariant measures (and, in particular, quasi-stationary 
distributions) of subcritical BGW processes, yielding for the first time a complete characterization 
of these measures involving an explicit formula. 
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